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We present a modified finite temperature Lanczos method for the evaluation of dynamical and 
static quantities of strongly correlated electron systems that complements the finite temperature 
method (FTLM) introduced by Jaklic and Prelovsek for low temperatures. Together they allow 
accurate calculations at any temperature with moderate effort. As an example we calculate the 
static spin correlation function and the regular part of the optical conductivity o re3 (u)) of the one 
dimensional Hubbard model at half-filling and show in detail the connection between the ground 
state and finite temperature method. By using Cluster Perturbation Theory (CPT), the finite 
temperature spectral function is extended to the infinite system, clearly exhibiting the effects of 
spin-charge separation. 

PACS numbers: 78.20.Bh, 71.10.Fd, 75.25,+z 



The Finite Temperature Lanczos Method (FTLM), in- 
troduced by Jaklic and Prelovsek Q , has in recent years 
allowed the precise calculation of thermodynamic quan- 
tities of strongly correlated systems. It has been applied 
to the t-J model for the cuprates |2|, |[ [| [?], |8| and 
vanadates || orbital t-J model p"lf, Kondo lattice 
model Hcisenberg model jll| and static properties 
of the Hubbard model jl4|, [l5j . In principle this method 
can be applied at all temperatures, but at low temper- 
atures the required number of random samples is very 
large. FTLM is restricted to small systems and particu- 
larly at low temperatures, finite size effects become im- 
portant. They can be overcome, at least for properties 
derived from the single particle Green's function, by us- 
ing Cluster Perturbation Theory (CPT) (T^, which 
provides infinite system results with remarkable accuracy 
p8| . However, CPT requires finite cluster Greens func- 
tions G a b{io) for all sites a,b, increasing the number of 
required Lanczos runs and matrix elements by a factor 
equal to the number of lattice sites. It therefore requires 
a sufficiently fast low temperature method. 

In this paper we present a modified finite temperature 
Lanczos method which allows to calculate properties for 
large Hilbert spaces at low temperatures that are not ac- 
cessible by the existing method. We show that a combi- 
nation of our low temperature Lanczos method (LTLM) 
with the FTLM allows an accurate calculation of thermo- 
dynamic properties at any temperature with moderate 
effort. 

Let us first present the existing FTLM. For the case of 
a static expectation value of an operator O 

N N 

(O) = -J2(n\Oe-? H \n), Z = ^(n| e -^|n), (1) 

n n 

with [3 = 1/T (fcs = % = 1) and a sum over a complete 
orthonormal basis set |n), the FTLM approximation is 
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with a sum over symmetry sectors s of dimension N s , 
and R random vectors \r) in each sector M is the 
number of Lanczos steps. For each random vector |?*), 
a Lanczos procedure is performed, yielding M eigenener- 

(r) . . (r)\ 

gieS Em and corresponding eigenvectors y&m ). The trace 
in (j|) requires N s states, while very accurate results can 
be obtained via (||) even for a drastically reduced number 
of Lanczos steps M N s and a partial random sampling 
of the Hilbert subspaces R <C N s . 

For dynamical correlation functions C(t) = (A(t)B), 
FTLM calculates 

1 AT R M 

s r i,j 

x (r\^)(^\A\¥;>)(^\B\r), (3) 
Here, an initial vector 

\^ ] )=B\r)/^{r\B^B\r). (4) 

is used to generate additional cigenencrgies ij and 

eigenvectors l^^') from that part of the Hilbert space 
onto which the operator B projects. Hence the term 
(J&^\B\r) in (|^) becomes sufficiently large. 

Let us now consider the behavior of ^ and (||) in the 
limit T — > 0. In this case only the ground state |^o) 
contributes and we get 

R R 

(O) = ^<*o|0|r)<r|* )/E(*ol r -)( r 'l*o) ( 5 ) 

r r 

and similarly for (||). Thus the ground state result will 
suffer from severe statistical fluctuations, although the 
exact (Lanczos) eigenvector |^o) is reached with every |r) 
and one random vector should be sufficient. Yet, FTLM 
gets worse with decreasing temperature T . 
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The modifications we present in this paper are designed 
to overcome this limitation. Let us introduce the method 
for a static expectation value (Q). We use a symmetric 
form 

1 N 

<0} = ^£He-^Oe-^V). (6) 

n 

As before, we approximate the trace by random sam- 
pling, but now we insert the approximate eigenbasis ob- 
tained by the Lanczos procedure twice, initially obtaining 

1 AT R M 

s r i,Z 

x (rl^X^IOI^X^Ir) (7) 

The partition function Z is calculated in the same way as 
in standard FTLM. The behavior in the limit T -> is 
now different. If only the ground state \ty ) contributes, 



(0) becomes 

(O) = 2<* |r)<r|*oX*o|0|*o)/X;<*o|r)<r|*o) 

r r 

= (#o|0|*o>. (8) 

In agreement with ground state Lanczos, one random 
vector suffices for the ground state expectation value. 

If we compute the numerator in (0) and Z sepa- 
rately, both suffer from pronounced statistical fluctua- 
tions, which however cancel exactly at T = as shown 
in (||). For finite T the fluctuations in numerator and 
denominator do not cancel exactly but they are still 
strongly correlated. Separate error analysis for both 
terms would overestimate the statistical noise. These cor- 
relations are taken into account by employing a Jackknife 
technique [po| . 

For dynamical correlation functions, a straight-forward 
variant of (0) suitable for low temperatures is 



C(t) = §£^££e-*«'^%^^ (9) 

s r ijl 



In order to span the relevant subspace of the Hilbert 
space, we now choose initial vectors |$q ) ^ ) f° r 

the second Lanczos run. With M such second Lanczos 
runs, the numerical effort would be much higher than for 
FTLM. For low temperatures, it can be reduced, since 
only the low lying states contribute to the expectation 
values. We consider only states below a cutoff energy 
E c , defined by 

e -ME c -E ) < ^ (1Q) 

where e c defines the accuracy of the approximation, (3 C is 
the minimal inverse temperature considered and the cal- 
culation will be accurate for all > f3 c . We thus proceed 
as follows: For each random start vector |r), we perform 
an initial Lanczos run with M iterations. For each of 
the M c states IV^) with energies below E c , we then cal- 
culate an initial vector |$q ) 00 ^l^i ) an d perform a 
second run with M Lanczos iterations, obtaining an ap- 
proximate eigenbasis |^ r '^). Using these basis sets, the 
final form of LTLM is the same as (0) and ([)]), with 

and replaced by J2i i an d > respectively. 

Memory requirements of our method are the same as 
for standard FTLM, but the CPU time requirements dif- 
fer significantly. CPU time is mainly determined by the 
number of matrix elements that have to be calculated. 
In the case of static expectation values these are M for 



FTLM and for LTLM for each random vector. There- 
fore both methods reach equivalent CPU time require- 
ments per random vector when M c m \pM. 

For dynamical correlation functions, the number of 
matrix elements to be calculated in the second Lanczos 
run is M 2 for FTLM and Af c 2 for LTLM. For LTLM we 
have to perform M c second Lanczos runs, but only one 
for FTLM. Thus we have similar CPU time requirements 
per random vector for both methods when M c M 3 . In 
the limit T^Owe have M c = 1, and for R = 1, LTLM 
is comparable to the ground state Lanczos technique. 

For both methods, CPU-time is proportional to R, the 
number of random vectors. But, by design, far less ran- 
dom vectors are needed for the LTLM than for the FTLM 
at low temperatures for a given accuracy. 

All considerations so far have been done without re- 
garding reorthogonalization of Lanczos vectors. This pro- 
cedure becomes important for M > 150, where numerical 
round off errors become significant, and significantly in- 
creases CPU requirements. 

Let us now demonstrate the method for the calculation 
of static and dynamical properties of the one dimensional 
Hubbard model, with hamiltonian 

H = ~ l X^lci+M +h.c.) + uJ2 n n n n- ( n ) 

Z,£T 1 

We specify energies in units of t. As an example we cal- 
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culate the static spin correlation function C\ — (S*S* +1 ) 
on a 12 site chain with periodic boundary conditions 
at half filling (n = 1). The number of basis states is 
N = 2 704 156. Symmetry sectors are specified by mo- 
mentum k and total spin S z . The largest sector S z = ±1 
has 52 272 basis states. The sector S z = is further 
reduced due to spin up/down symmetry. 
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FIG. 1: Spin correlation function Ci = (S^ Sf+i) for the ID 
Hubbard model on a 12 site chain with periodic boundary 
conditions at U=8t and n=l. Solid: LTLM with (3 C = 1, 
e c = 0.01, M = 100, and R = 25. Dashed: Two independent 
runs of FTLM with M = 100 and R = 25. Inset: Deviation 
of Ci in the high temperature region beyond (3 C . Here R = 50 
in both cases. 
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FIG. 2: Relative statistical errors of LTLM (solid) and 
FTLM (dashed) with R = 25 in both cases. Other parameters 
as in fig. [j]. The error of LTLM at T = O.Olt is A rei = 10" 8 . 

In figs. | and § the convergence and statistical errors 
of LTLM and FTLM are compared at equal computa- 
tional effort, with R = 25 random samples per symmetry 
sector each, corresponding to sampling of S- w 0.05% 
of the largest Hilbert subspace. At low temperatures, 
our method provides results which are orders of magni- 
tude more precise than from standard FTLM, and which 
connect smoothly to the ground state properties. We 
checked that for larger R, there is no systematic drift for 
either method, and the FTLM results converge towards 
those of LTLM. At intermediate temperatures, the sta- 
tistical errors of LTLM increase and become similar to 
those of FTLM. Finally, considerably beyond the chosen 
cutoff-temperature 1//3 C , LTLM is no longer applicable, 
and begins to show a systematic deviation. 

Both FTLM and LTLM provide results for a range 
of temperatures from a single calculation. For FTLM 




FIG. 3: Regular part of the optical conductivity of the Id 
Hubbard model on a 12 site chain with periodic boundary 
conditions at U=6t and n=l. Left panel: LTLM calculations 
with p c = 3, £ c = 0.01, and R = 40. Right panel: FTLM 
calculations with R = 50. Number of Lanczos steps M — 100 
and additional broadening of a — 0.1. Dots mark the zero 
line. Only selected errorbars are shown. For curves without 
errorbars, the errors are smaller than the line width. 



this range is limited towards low temperatures by sta- 
tistical errors. For LTLM, it is limited by the chosen 
cutoff-temperature 1//3 C . Therefore a combination of 
both methods provides precise results for all tempera- 
tures with moderate effort. 

As an example of dynamical correlation functions we 
calculate the regular part of the optical conductivity, 
given by the current-current correlation function 
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dte iut {j(t)j), (12) 



with the current operator j = it JT -( c io- c »+L°' — h- c -)- ^ n 
fig. |^, we show results with approximately the same CPU 
time for both methods. Slightly above the ground state, 
at (3 = 40, LTLM approaches the exact ground state re- 
sult ||l], |2^]. For intermediate temperatures [3 — 10, 5, 3, 
slight statistical fluctuations occur. By comparison to 
FTLM we see that (3 — 1 < f3 c is indeed beyond the va- 
lidity of this calculation. We also checked the accuracy 
of the results by using M — 200 Lanczos steps instead of 
AI = 100 yielding the same LTLM spectra within statis- 
tical errors. 

In contrast, FTLM suffers from strong statistical fluc- 
tuations at small temperatures. Errorbars are very large 
and regions occur where a re9 (uj) becomes negative, a 
clear indicator that we did not use enough random vec- 
tors for FTLM. As expected from our consideration of 
static expectation values, errorbars of FTLM get smaller 
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for higher temperatures. As for LTLM we did calcula- 
tions with M = 200, yielding the same curves within 
errorbars but leading to a better convergence at the high 
frequency side of the spectrum. 





FIG. 4: Spectral function A(k,w) obtained by LTLM for the 
ID Hubbard model at U = 8t, n = 1, and f3 = 10. Parame- 
ters: R — 30, p c = 3, £ c = 0.01. Upper panel: 8 site chain 
with periodic boundary conditions. Lower panel: CPT result 
based on 8 site clusters. 

As mentioned in the beginning, at low temperature 
finite size effects become important. At least for proper- 
ties derivable from the single particle Green's function, 
these effects can be overcome by using CPT 16, p~7| . 
In fig. ^ we show the finite temperature spectral func- 
tion A(k, <J) obtained on a finite size lattice with peri- 
odic boundary conditions, and the infinite lattice result 
obtained by CPT, which makes use of all Greens func- 
tions on the finite lattice as calculated by LTLM. In the 
latter a smooth structure consisting of several branches 
can clearly be seen with spin-charge separation at k = 
visible in the lower part of the spectrum jl(| |2j| . On the 
finite size cluster, however, this structure is not evident 
as it exhibits more discrete excitations. Further work on 
finite temperature CPT is in progress pl| . 

In conclusion, the method presented in this paper gives 
an accurate connection of the exact ground state Lanc- 
zos method and the established FTLM. Using LTLM at 



low and FTLM at higher temperatures makes it possible 
to calculate static and dynamical properties of strongly 
correlated systems from T = up to T = oo with very 
good accuracy and rather small numerical effort. 
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